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Abstract 

In this work, we continue the modelling of the digestion in the small intestine, 
started in a previous article, by investigating the effects of dietary fibre. We recall 
that this model aims at taking into account the three main phenomena of the diges- 
tion, namely the transit of the bolus, the degradation of feedstuffs and the absorption 
through the intestinal wall. In order to study the role of dietary fibre on digestion, 
we model their two principal physiochemical characteristics which interact with the 
function of the small intestine, i.e. viscosity and water holding capacity. This leads 
us to consider some features of digestion which have not been taken into account 
previously, in particular the interrelationship between the evolution of dry matter and 
water in the bolus. The numerical results are in agreement with the positive effect 
of insoluble dietary fibre on the velocity of bolus along the small intestine and on 
its degradation. These results highlight the negative effect of soluble dietary fibre on 
digestion. Although, this model is generic and contains a large number of parame- 
ters, to our knowledge, it is among the first qualitative dynamical modelling of fibre 
influence on intestinal digestion. 

Key-words : Modelling of intestinal digestion, water holding capacity, soluble dietary fibre, in- 
soluble dietary fibre. 

1 Introduction 

Digestion in the small intestine can be described through three main phenomena : transit of the 
bolus along the small intestine, degradation of macromolecules into smaller ones and absorption 
through intestinal wall. Taking into account these phenomena, the authors have presented in [1] 
a generic model of digestion in which the bolus include only one category of macromolecules 
(carbohydrates, proteins or lipids) and water. 

However mixing these nutrients influences the digestion process through interactions between 
molecules. In order to improve this model and to make it more realistic, we should consider the 
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effects of such interactions, and we have decided to do so by first including dietary fibre in the bolus 
because of their significant role on the digestion. One of the key properties of fibre is its water 
holding capacity and this leads us to investigate the role of water in the digestion processes. To do 
so, we distinguish dry matter and water in each substrate, we model water kinetic in correlation 
with the dry matter one and we take into account the impact of water in all the aspects of digestion. 

To be more precise on our approach, we first describe the main underlying assumptions which 
guide our modelling, then we show how they are translated into equations and finally numerical 
tests are performed for examining the effects of several hypothesis. In addition, we point out 
known or assumed mechanisms relating dietary fibre and digestion process. 

This paper is organized as follows : physiochemical characteristics of soluble and insoluble di- 
etary fibre are introduced in Section 2. In Section 3, main assumptions of the model are presented. 
Section 4 is devoted to the description of the composition of bolus, the chemical transformations 
of macromolecules and all the notations. Digestion in the presence of dietary fibre includes the 
modification of transport as well as degradation and absorption equations as described in Section 5. 
Section 6 is a comparison of the numerical results considered as in silico experiments. Finally, in 
Section 7, the model is discussed and perspectives are proposed. 
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INRA Center of Nouzilly and the Laboratoire de Mathematiques et Physique Theorique was ini- 
tiated within and supported by the CaSciModOT program (CAlcul SCIentifique et MODelisation 
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collaboration also takes place in a CNRS-INRA PEPS program "Comprehension et Modelisation 
du devenir de 1' aliment dans le tube digestif". This work is part of the PhD thesis of Masoomeh 
Taghipoor, financed by CNRS and INRA. 

2 Biological Background on Water and Dietary Fibre 

Digestion modelling requires the knowledges of the physiochemical properties of macro- 
molecules concerned by this phenomenon as well as mechanical and biochemical reactions ob- 
served for their degradation. 

Dietary fibre (DF) is usually defined as the sum of plant non-starch polysaccharides and lignin 
that are not hydrolysed by the enzymes secreted by the non-ruminant digestive system, but that 
can be partially digested by microflora in the gut. A main effect of fibre is to regulate intestinal 
degradation and absorption of nutrients as well as their transit along the gut. Physiochemical 
characteristics of fibre include viscosity, hydration, fermentability (mostly in the large intestine), 
adsorption or entrapment of nutrients and bulking effect. Each of these characteristics affects 
meaningfully the function of the gastrointestinal tract [2, 3] . These characteristics depend on the 
polysaccharides chemistry. One way to classify dietary fibre is based on their water solubility. 
Insoluble dietary fibre include cellulose, some hemicelluloses and lignin. The other is soluble 
dietary fibre such as viscous fibre which includes beta-glucans, pectins, gums, mucilages and 
some hemicelluloses [4, 5]. 

For monogastrics, most available nutrients are degraded and absorbed in the small intestine. 
At the beginning of duodenum the bolus consists of partially degraded feedstuffs and water. Once 
in the small intestine, mechanical and chemical digestion of feedstuffs make the nutrients avail- 
able to the organism. Enzymatic hydrolysis is the most important chemical reaction in digestion, 
which takes place in aqueous solution. Enough water is required for an efficient digestion process 
eventhough water/nutrient ratio are not precisely known. Furthermore, classification of dietary 
fibre through their water solubility and the impact of Water Holding Capacity (WHC) of DF on 
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digestion reveals the key -role of water on digestion. WHC is defined by B. Shneeman [3] as the 
ability of fibre source to swell when mixed with water and to hold water within its matrix. 

Soluble dietary fibre 

Soluble DF are believed to impact significantly digestion and absorption as well as bolus 
transport in the small intestine. The main physiochemical properties of soluble DF are viscosity, 
water holding capacity (WHC) and organic compound entrapment [6]. Soluble DF, because of its 
high viscosity, is generally associated with slow transit through the stomach and increasing of the 
small intestinal transit time [7]. 

Insoluble dietary fibre 

Insoluble DF acts primarily in the large intestine where, due to its WHC, increases faecal 
bulk, dilutes colonic contents and decreases mouth-to-anus transit time [6]. However, its effects 
on digestion and transit in the small intestine can not be neglected since insoluble DF affects the 
transit time in the small intestine through its laxative property. Recent studies have shown that the 
inclusion of a moderate level of dietary fibre improves the digestibility in chicks [8]. Therefore 
to obtain an optimal efficiency in nutrient utilization, Burhalter et al. [9] proposed to increase the 
ratio of insoluble to soluble DF. Moreover, the use of insoluble fibre in commercial broiler chicks 
improves the intestine morphological parameters and result in a better performance assumed to be 
connected to more efficient digestion and absorption processes [10]. Two hypothesis are proposed 
in order to study the influence of insoluble DF on nutrients digestibility in the small intestine : 
(i) insoluble DF increases the retention time in the stomach changing the nutrient profile of the 
bolus at the entry of the small intestine which could lead to a higher digestion and absorption, (ii) 
Physical characteristics of insoluble DF change the digestion process mostly through their capacity 
of swelling water and nutrients. These both hypothesis are either tested in the in silico experiments 
(cf. Section 6.1) or included in the equations (cf. Section 5). 

Water 

To have a better understanding of the role of dietary fibre, it is therefore required to study more 
precisely the evolution of water during digestion in the small intestine. For example, depending 
on the bolus composition, water absorption through intestinal wall or its secretion into the small 
intestine lumen could be observed. However, the evolution of water amount in the small intes- 
tine depends also on other components' kinetics within the bolus as described by (Reaction 1), 
(Reaction 2) and (Reaction 3) in Section 4. 

3 Key Model Assumptions 

In this section, the key assumptions for the model are presented. 

HI: Each component of the bolus (macromolecules, partially degraded macromolecules, nutri- 
ents and fibre) is represented mathematically as a portion of dry matter and a characteristic 
proportion of water. 

For example, "starch in a bolus" includes both dry starch and water used to maintain starch 
molecules in aqueous solution. The same is observed for the "disaccharides in a bolus" and "glu- 
cose in a bolus" combining smaller molecules resulting from starch hydrolysis associated with a 
specific level of water. In other words, a component C of bolus is represented as C dm + Wc 
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where C dm denotes the dry matter of C and Wc is the necessary amount of water to maintain it 
in a solution state. Moreover, we have assumed that the mass of Wc is proportional to the mass 
of C dm , i.e. equal to cC dm for some characteristic number c > whic represents the necessary 
amount of water proportion to maintain C in solubilized phase. Despite the presence of water in 
the bolus, a little amount of non-solubilized dry matter may be included in bolus, which is (of 
course) associated with c = 0. 

H2: Without DF, the bolus contains a single macromolecule and water. It is represented by a 
homogeneous cylinder with the constant length I. Including insoluble DF transforms this 
homogeneous bolus into an heterogeneous one by modifying the concentrations offeedstuffs 
and nutrients. 

A bolus in the small intestine is a viscous solution of dry matter and water: we assume that 
its volume is very close to the volume occupied by the water in the bolus. In other words, our 
assumption is that DM does not fill any volume (e.g. : solubilized sugar + water does fill the same 
volume as the water alone). Once the volume is known, [H2] allows to compute the radius of the 
bolus r(t), which impacts its movement along the small intestine. For a bolus with constant mass 
at the beginning of duodenum, [H2] influences the degradation by increasing the concentration of 
nutrients in the homogeneous part of bolus. Indeed the space occupied by insoluble dietary fibre 
is unavailable for the other macromolecules. 

H3: Digestion in the small intestine is due to volumic and surfacic transformations. Volumic 
degradation is the enzymatic hydrolysis of bolus components by pancreatic and exogenous 
enzymes inside the bolus while the surfacic one is the degradation by brush border enzymes 
on the internal wall of small intestine. 

Some additional facts should be pointed out : water facilitates the contact of the macromolecules 
with the brush border enzymes, enhancing the surfacic degradation. Increasing the water to dry 
matter ratio dilutes the bolus and decreases the volumic degradation. Both of these reactions are 
proportional to the mass ratio of concerned substrates. 

H4: The bolus movement along the small intestine is due to peristaltic waves. The efficiency 
of these waves are proportional to the radius of bolus, and inversely proportional to the 
distance from the pylorus. 

We use, in this article, the model of transport with an averaged velocity function as presented in 
[1 1]. The movement of the bolus in the model by this equation is due to a homogenized accelera- 
tion caused by the average effect of the peristaltic waves. Of course, there would no difficulty to 
come back to the model which takes into account all the frequent peristaltic waves. 

H5: The water in the bolus which is not hold by the macromolecules and DF through WHC i.e. 
"available water", decreases the viscosity of the bolus and facilitates its movement. Due to 
osmotic type equilibrium, the concentration of this "available water" tends to reach a fixed 
ratio. 

In other words, "Available water" reduces the friction caused by the bolus contact with the intesti- 
nal wall. 

H6: Dietary fibre modifies the bolus evolution through its WHC by holding the water in its 
matrix, and therefore modifying the volume of the bolus. Soluble DF decreases the efficiency 
of peristaltic waves. 

By their WHC, dietary fibre holds a significant quantity of water in the bolus and therefore keep 
the volume and the radius of the bolus higher. It is worth pointing out that soluble dietary fibre 
change the consistence of the bolus in the sense of making it more jelly, implying the decreases of 
efficiency of peristaltic waves. 
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4 Physiological Aspects and Bolus Composition 



Different steps of mechanical and chemical transformations are detailed below. We present 
also the composition of the bolus and their interactions. 

• Non degradable substrate A na - : The quantity A na - represents the mass of macromolecules 
which is not degradable by endogenous enzymes of the digestive tract. 

• Non solubilized substrate A ns : The quantity A ns represents the mass of macromolecules 
which is not accessible to enzymatic hydrolysis. In presence of a sufficient quantity of 
water, A ns is transformed into A s . 

• Solubilized substrate A d s m : The quantity A d s m is the mass of dry substrate in solution 
state. It is called solubilized substrate and it is assumed that one unit of A dm requires W s 
units of water to remain solubilized. Recalling [HI], W s represents the required mass of 
water to solubilize Af 71 . This quantity depends on the properties of each macromolecules. 
The mass of W s is assumed to be equal to the mass of aA d s m where a represents the ratio 
of water associated with A dm . For simplification purposes, the mix of A dm and water is 
represented by A s . 

• Intermediate substrate Bf^ : The quantity Bf^ is the mass of dry intermediate substrate 
obtained from the degradation of A s by volumic transformation [H3]. It is solubilized and 
Wi n t represents the required amount of water to maintain solubilization. The mass of Wi n t 
is assumed to be equal to the mass of fiBf^ where (3 represents the ratio of water associated 
with Bf^. For Bi nt = Bf^ + W{ n t, volumic transformation is represented as 

A s + enzymes ->■ B int + (W s - W int ) (Reaction 1) 

Depending on the value of W s and Wi n t, the amount of W s — Wi n t of water can be released 
or hold in the bolus. 

• Absorbable nutrients B d ^ : The quantity B d ^ is the mass of dry absorbable nutrients 
obtained from surfacic reactions (cf. [H3]). For B abs = B d ^ + W a6s , the surfacic transfor- 
mation is defined as 

A s + enzymes ->■ B abs + (W s - W a bs) (Reaction 2) 

B int + enzymes ->• B ahs + (W int - W abs ) (Reaction 3) 

W a bs is the required amount of water to maintain solubilization and its mass is assumed to 
be equal to 7 B^™ where 7 represents the ratio of water associated with the B%£. 

• Soluble and insoluble dietary fibre and Ff™ ol represent the soluble and non-soluble 
dry dietary fibre respectively. The main property of dietary fibre presented in Section 1 and 
hypothesis [H5] is its water holding capacity. 

F sol = W sol + F^ 

Finsol — ^^insol -^insoi 

where the mass of W so i {Wi nso i) is assumed to be equal to \ s Ff^ {\F^ ol ) for \ s and 
\ which represent the ratio of water associated with Ffj^ and F^ ol respectively . As 
described in Section 1 , DF is not degradable by endogenous enzymes of the small intestine. 

The following diagram shows the different transformations inside the bolus. 
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Figure 1: Physical and chemical transformations inside the bolus and included in the model are 
represented in this scheme. "B.B" stands for brush border enzymes and "per" stands for pancreatic 
ones. 

• Dry Matter : Total amount of Dry Matter substrate in the bolus is therefore 

DM = A nd + A ns + Ai m + B% + BiZ + F% + F?™ 



insol ■ 



• Water : Impact of dietary fibre on digestion is closely linked to their WHC capacity. 
Though, water evolution in the bolus has to be described to understand effects of DF on 
digestion. 

Total water Wtot in the bolus comes from three main sources 

(i) Wf ee d : water incorporated naturally in feedstuffs (e.g. : one gram of wheat contains 
12 % of water). The amount of Wf ee d coming out of the stomach is assumed to be 
proportional to the ingested dry matter DM, Wf ee d = Kf ee d ■ DM. 

(ii) W sec : Water included in the endogenous secretions of saliva and stomach which is 
also assumed to be proportional to the ingested dry matter, W sec = K sec ■ DM. 

(Hi) WdHnk '■ Drunk water is assumed to be independent to the quantity of DM. 

Total water included in the bolus is therefore defined as the sum of Wf ee d, W sec and WdHnk 

Wtot = Wf ee d + W sec + W drink- 

Thereby, "available water", W, as presented in [H6], is defined as the difference between 
Wtot and the water associated with DM for maintaining the solution state. In term of mass, 
the quantity of "available water" in the bolus at each time is 

W = Wtot ~W S - W int - W abs - Wsol - W inso l . 



• Mass of bolus M : The total mass of bolus, M, is given by M = DM + Wtot- 
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• Volume of bolus V : To define the bolus volume, as explained in [H2], the volume of each 
substrate in solution is assumed to be the same as the volume filled by water associated with 
that substrate, i.e. the volume of A s is equal to W s /p w , where p w = 1 is the water density. 
The volume of bolus is therefore represented as 

V = W tot = 7rr 2 e Pw = irr 2 £. 

Since the length of the bolus is assumed to be fixed, the volume evolution leads to compute 
the radius r(t) of bolus at each time. Consequently its surface is written as S = 2nr£. 

In the following section, the different properties of dietary fibre on the model of digestion 
are taken into account. 



5 Model Equations 



To include WHC property of soluble DF in the digestion model, the mix of F^™ and W so i is 
assumed to form a viscous gel in the bolus. Therefore in our model the mass of "available water" 
in the bolus is reduced. Moreover, viscous fibre enhances motility but decreases transit rate, since 
it resists propulsive contractions [12]. This resistance to peristaltic waves is described through a 
new notion called efficient radius of bolus called r so i. As described in [H2], the volume filled by 
soluble DF is W so i = \ s F so i, then 

Definition 5.1. The efficient radius of bolus is defined as 

r sol = ^{W tot - W so1 )/2kL 
In the same way, the efficient surface of bolus is described as S so i = 2nr so iL 

This definition is used to describe the decrease of the surfacic degradation and absorption 
caused by soluble DF in the model. 

One of the hypothesis in the first model of digestion in [1] is the bolus homogeneity. The 
mass concentration of each component of the bolus is assumed to be its mass divided by the total 
mass M of the bolus. To model the digestion in presence of insoluble DF, new notions are defined 
because of heterogeneity of bolus as described in [H2]. The volume filled by insoluble DF (the 
mix of Ff™ ol and Wi nso i), is assumed to be unavailable to the macromolecules of feedstuffs in the 
bolus. In the digestion model, this hypothesis is taken into account by the following definition. 

Definition 5.2. The apparent concentration of different substrates in the bolus is represented as 

A dm r>dm p>dm 

r Admi s rndmi _ int [n£im| _ abs 

1 s ' ~ l\/T TP ' y int \ ~ jiyf Tp i i^absl ~ n,r t? 

- r inso i M - £ inso l M - t inso x 

A ns .„_. W .. . A nd 



[^]= Z , [W] = - " - , [A 



M-F insol ' M-F insol ' ™ l M-F msol 

The degradation of macromolecules A s and Bi nt as well as the absorption of the nutrients 
B a f, s are affected by WHC property of the insoluble DF through this definition. Figure 2 shows the 
regions of bolus which are filled by insoluble DF and therefore unreachable by the macromolecules 
and nutrients. 
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Figure 2: The distribution of insoluble fibre in the bolus as assumed in the model. The apparent 
volume V app (see Definition 5.3) is the white part of the cylinder. 

Integrating the insoluble DF in the bolus changes also the region reachable to volumic degra- 
dation (see Figure 2). According to hypothesis [H2], Wi nso i is the volume filled by insoluble 
fibre. 

Definition 5.3. The apparent volume of bolus V app is therefore defined as 

V app = Wtot — Winsoi. 

The above considerations are taken into account in the following steps of digestion described 
below. 



5.1 Transport of bolus 

The averaged equation of transport of bolus introduced in [11] reads 

d 2 x _ _ 1 dx cq + ci r K visco dx 

dt 2 ~ T{ C dt> a + bx [W] dt ' 

where r is the mean effect of the pulses by unit of time, c is mean velocity of peristaltic waves, x 
the position in the small intestine. Taking into account the properties of dietary fibre, this equation 
changes to 

d 2 x _ _! dx cq + cir soi K visco dx 

dt 2 ~ T[ C dt' a + bx [W] df { ' 

The bolus movement described by this equation depends on its position in the small intestine 
and on its efficient radius. Moreover, the acceleration is slowed down by a viscosity term which 
depends on the available water. 



5.2 Volumic transformation 

Volumic transformation presented by (Reaction 1) in Section 4 is the degradation of Af™ due 
to pancreatic and exogenous enzymes resulting in production of Bf%£. Evolution of A^ m by this 
transformation is represented by 

dA dm d 
^ = — k vo i(x){A s ]V app 

where k vo i(x) takes into account the enzymatic activity which is a function of bolus position 
at each time t. The term [y^ m ]V app describes the dependence of volumic degradation on the 
concentration of A^ m at each unit of apparent volume i.e. the volume filled by the insoluble DF is 
not accessible to the enzymes and macromolecules. 
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Consequently, integrating the insoluble fibre in the bolus influences the volumic transformation 
by increasing the substrates concentration via the Definition 5.2 and by changing the volume and 
using V app introduced in Definition 5.3. 

As described in Section 5.1, integrating the soluble fibre modifies the velocity of bolus along 
the small intestine and therefore the distance travelled at each time x(t). Consequently, it influ- 
ences the volumic degradation through the function k vo i{x(t)). 

This reaction takes place in a solution and each unit of A dm degraded by this equation causes 
the release of "available water" W s = aA dm \ 

The volumic production of intermediate substrate Bf^ is the result of degradation of A dm 

dF> dm ^ 

J^ 1 = kvoi( x )[A s m ]V ap p. 

Each unit of produced Bf^ requires and uses Wi n t = PBf^ to maintain the solution state. Ac- 
cording to (Reaction 1), the result of volumic transformation is the consumption or release of 
"available water" W. Thereby, volumic evolution of water is represented as 

d ^ = ~k vol ( x )[Af n ]V app , (2) 

where k voi = (a - /3)k vo i. 

Soluble fibre can be hydrolysed by exogenous enzymes e exo (cellulase, hemicellulases, ...) 
ingested by food 

Fgol + e exo -> B int + {W Q i - W int ). (Reaction 4) 

Hydrolysis of F so i by exogenous enzymes follows the same evolution as the volumic transfor- 
mation of A dm 

jrpdm 

= -k s e exo ph(x)[F^]V app 

where [F^™] = F so i/(M — Fi nso i) and ph(x) is the exogenous enzyme activity along the small 
intestine. This reaction produces the intermediate substrate Bf^ 

j r>dm 

—^ = k s e exo ph(x)[F^}V app . 
The amount (W so i — Wi n t) is released by (Reaction 4) and modifies the evolution of water 

^ = ... + k s e eX oph(x) [F^] V app (3) 

where k s = (X s — f3)k s . 

5.3 Surfacic transformation 

Surfacic degradation is the last step of transformation of macromolecules in the small intestine. 
The produced nutrients by this degradation are then absorbed through intestinal wall. Surfacic 
degradation depends on the fraction of A s on the surface of the bolus represented by [A^J^o/ 
therefore 

fjAdm 

= -k surf [A s ][W]S sol 

where k sur j stands for the rate of surfacic degradation of A s and the efficient surface S so i has been 
defined by Definition 5.1. Moreover, it is assumed that the brush-border enzymes are always in 
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excess in the small intestine. Surfacic degradation of Bf^ follows the same process as for Af 11 . 
Therefore for [Bi nt ] defined by Definition 5.2, we have 

Jj* = -k SU rf[Bint\[W]S so i. 

where k sur f stands for the rate of surfacic degradation of Bi nt . Evolution of water in the bolus 
is influenced by the surfacic degradation, i.e. the quantity of water consumed (or released) by 
(Reaction 2). Therefore the surfacic evolution of water is 



dW 
~dT 



- + - l)kurf[Bint] + (a - i)k surf [A s ^j [W]S sol . (4) 



5.4 The equilibrium between A s and A 



ns 



Modifications of feedstuffs in the stomach by the enzymes and water change most of the A ns 
into A s = Ag" 1 + W s and makes them accessible to intestinal enzymes. 

However, for some feedstuffs, the bolus may contain A ns at the beginning of small intestine. 
In this case, the digestion of bolus contains also the transformation of A ns into A s . Mixing with 
bile acid for lipids and producing the micelles, denaturing for the proteins and adding water and 
solubilization for the dry starch are examples of the transformation of A ns into A s in the small 
intestine. 

The solubilization of A ns which takes place in the presence of enough quantity of W and 
results in the production of A s , is a phenomenon taken into account in the model. Solubilization 
is reversible and lack of water may cause production of A ns releasing W in the bolus. 

Thereby, the balance is assumed to be reached when 

A s = fi([W])A ns 

for n which is an increasing function of [W] . If k equ i stands for the rate of turning back to equi- 
librium then the dynamical equilibrium may be defined 

A A 

-^ = -k equi {fi([W])A ns -A s ) (5) 

and therefore 

= k equi (n([W])A ns - A a ). 
The variation of water quantity caused by the equilibrium may be represented as 

^ = ... + ak egui ((i([W])A ns - A s ). (6) 

5.5 Pancreatic and biliary secretions 

Pancreatic and biliary secretions consist of a solution of nutrients and enzymes which do not 
contain available water W. In fact, water included in this solution is assumed to be associated 
with nutrients and enzymes to keep them solubilized. Modelling details on these secretions could 
be seen in [1]. Adding dietary fibre increases the quantity of pancreatic secretions. However this 
point is not yet included in the model. 
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5.6 Absorption through intestinal wall 

Absorption of nutrients through intestinal wall depends on their concentration on the inner 
surface of intestinal wall 

dB dm 

abs U \ T>dm] n 

7, — ••• ^abs[- D abs\ D sol 

The passage of nutrient through intestinal wall releases the associated water, thus 

d ^ = ... + 1 k abs [B±]S sol (7) 
where k a b s represents the rate of absorption through intestinal wall. 

5.7 Water equilibrium 

Water equilibrium was already taken into account in [1]. The assumption was that [W] tends 
to reach a fixed ratio (10%), suggesting the equation 

d M = -k w ([W]-0.1) (8) 

where [W] = W(t)/M(t), M(t) representing the bolus mass. The superposition of the equations 
(2), (3), (4), (6), (7) and (8) provides the equation describing the evolution of W along the small 
intestine. 

The evolution of bolus mass is represented as 

dM _ (n , u dA t m , ,n , u dB?% , u dB± , u dF^ dA ns dW 

each term of the above equation is replaced by its expression, therefore we obtain 

W = whv { ~ kw{W " °- 1M) " k ^[BiZ\S S oi). (9) 

The variation of bolus volume depends on the absorption or secretion of "available water" and 
endogenous secretions in the small intestine i.e. 

therefore 

— = -k w (W - 0.1M) + Secretions. 



6 Numerical Simulations 

A thorough examination of the effects of the different parameters of the model on transport, 
degradation and absorption is carried out by Scilab software. In these in silico experiences, the 
presence of dietary fibre F so i and Fi nso i in the bolus, the variation of the initial values of bolus and 
the sensitivity of the model on the parameters a, (3 and 7 (cf. Section 4) are investigated. 

The mass of the bolus at the entry of the small intestine is assumed to be fixed in all our fol- 
lowing experiences. The bolus at the entry of the small intestine contains A n &, A ns , A s , Bi nt , B a b s 
and W. When studying the influence of dietary fibre on digestion, the non degradable substrate 
A nd is replaced by F so i or F inso i. 
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6.1 Influence of dietary fibre on intestinal absorption 



As described in Section 1 , the positive effect of insoluble DF in digestion may be due to two 
main reasons : (i) modification of the composition of bolus due to the increase of retention time 
in the stomach, (ii) modification of the bolus physical characteristics. 

The effects of these two cases on the digestion model are studied in this section. 

Influence of the modification of the bolus in the stomach 

Including insoluble DF in the bolus delays gastric emptying. The direct effect of this phe- 
nomenon is to increase the solubilization in the bolus and to start partially the degradation. We 
studied the effects of such a change in the initial conditions for our model. To this aim, two nu- 
merical experiences are carried out : (a) the increase in the ratio of A s to A ns and (b) the increase 
in the ratio of Bi nt to A s when A ns = 0. 

(a) Our first experience consists in increasing the ratio of A s to A ns in the bolus at the entry 
of the small intestine. The value of absorbed dry nutrients at the end of the small intestine 
does not vary meaningfully. Table 1 shows the absorbed dry nutrients to DM ratio at the 
end of ileum x = 17 for different for the different A^ m to DM ratio at the beginning of the 
small intestine. Despite the variations in the ratio of A s to A ns , the equilibrium between A s 





x = 


= 




End of the small intestine x — 17 


A^.DM- 1 < 


;%) 


A^.DM- 1 ( 


%) 


Absorbed dry nutrients to DM ratio (%) 







85 




56 


42 




42 




57 


85 









58 



Table 1 : The proportion of absorbed nutrients to DM for different scenarios of solubilization in 
the bolus at the entry of duodenum. 



and A ns defined by Equation (5) is reached quickly (see Figure 3). The sensitivity analysis 
shows also that the value of absorbed dry nutrients is almost independent to the variations in 
the ratio of A s to A ns . However, this results depend on the choice of equilibrium rate k equ i, 
a small value of k equ i may decrease the difference between the final values of absorbed 
nutrients. 

(b) In the second experience, the modification in the stomach is assumed to result in the trans- 
formation of all A ns in the bolus into A s and additionally the production of Bi nt . Variations 
in the ratio of A s to Bi nt inside the bolus at the entry of the small intestine are tested. Table 
2 shows the variation of absorbed dry nutrients at the end of the small intestine as a function 
of the initial value of Bi nt . 

Numerical results shows the increase in absorbed dry nutrients when the ratio of Bi nt to A s 
is increased. 

Direct effect of DF on the function of the small intestine 

Besides the modification of the bolus in the stomach, presence of insoluble DF changes also 
the physiochemical characteristics of bolus (Section 5). 

To observe the effect of DF in the model of digestion, value of insoluble and soluble DF was 
increased from 1 g to 5 g in a bolus of 120 g. 
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Time (hour) 

(c) A* m = 0, A ns = 36 



Figure 3 : The equilibrium Ag m — A ns is reached quickly for different initial value of A. 
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Table 2: The relation between the absorbed dry nutrients at the end of digestion and the different 
initial values of BflJ}. 



Figure 4 shows that the presence of insoluble DF promotes intestinal absorption, however 
this increase in absorbed dry nutrients is not meaningful. The results in this figure, show that the 
increasing of the value of soluble DF decreases the quantity of absorbed dry nutrients and increases 
the final total mass. 

Therefore, numerical simulation shows that the positive effect of insoluble DF on the amount 
of absorbed dry nutrients is mainly due to the modification of the bolus in the stomach. However, 
the effect of the interaction of bolus with the small intestine as described in Figure 4 may not be 
neglected. 

Time of intestinal transit in the presence of DF 

Numerical results of transit time in presence of DF are presente in Table 3. 
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(a) Absorption versus F,- ns (b) Absorption versus F s 

Figure 4: Change in the final mass of bolus and absorbed dry nutrients for different amount of 
DF in the bolus at the entry of duodenum. 



These results show that integrating insoluble DF in the bolus decreases the time of intestinal 
transit from 5 h for a bolus of 120 g without insoluble DF to 3, 9 h for a bolus of the same mass 
which contains 5 g of DF These results are consistent with published values. The experiences 
done by Wilfart et al. [13] have shown that increasing dietary fibre content reduced or tended to 
reduce the mean retention time in the small intestine. 

These results show that integrating soluble DF in the bolus increases the intestinal digestion 
time from 5 h to 6, 7 h illustrating the effect of viscosity due to soluble DF on transit time. 
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Table 3: Intestinal transit time for the different quantities of F 
entry of duodenum 



F,„i or Fi 



insol 



in the bolus at the 



6.2 Water associated to dry matter 

To study the influence of the quantity of associated water on digestion, variation of the values 
of a, (3 and 7 have been tested. We observe their influence on digestion and specifically on the 
absorbed dry nutrients and A s -A ns equilibrium. Two main simulations are carried out : an uniform 
water content for A s , Bi n t and B a b s i.e. a = (3 = 7 and a a non uniform one. 

The objective of these experiences is to understand how the different values of a, (3 and 7 
influence the digestion in our modelling. 

6.2.1 Uniform water content for A s , B int and B abs 

The value of a = (3 = 7 varied from 1 to 4 in the model presented in section 5. Our objective 
is to observe its effect on the value of absorbed dry nutrients as well as on the final mass of bolus. 
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Numerical results presented in Figure 5 show the negative effect of this increase on the ab- 
sorbed dry nutrients. Increasing the quantity of water (a, /3 and 7) associated with the dry feed- 
stuffs (Ag m , B^ t and B%™) in our model, dilutes the bolus and decreases the volumic degradation, 
it decreases also the quantity of dry nutrients in contact with the internal surface of the bolus. These 




Figure 5 : Dependence of the absorbed dry nutrients and the final mass of bolus at the end of the 
small intestine on the value of a, f3 and 7. 

results seems to be consistent with the reality, in fact, the more water is presented in the bolus, the 
less (pancreatic and brush border) enzymes and molecules are in contact. 

We are also interested by the effect of these variations on the A^ m — A ns equilibrium defined 
by Equation (5). The results are shown in Figure 6. The equilibrium is almost reached in the four 
experiences, however the choice of k equ i can change the necessary time to reach the equilibrium. 
A significant production of A ns is observed in Figure 6.d because of the lack of the available water 
at the beginning of the small intestine. 

6.2.2 Non-uniform water content for A s , B int and B abs 

In the second experience, the quantity of absorbed dry nutrients at the end of the small intestine 
and the numerical results of A s -A ris equilibrium for different values of a, (3 and 7 were observed. 
The choice of the values of A s , Bi nt and B a b s is based on the hypothesis that the value of (3 is 
always between the values of a and 7 because of the molecule size of Bf^. Even if the longest 





Absorbed dry nutrients to 
DM ratio (%) 


Retention time in the small 
intestine (h) 


a = 1,0 = 2, 7 = 3 


45 


5,2 


a = 2(3 = 2, 7 = 2 


50 


4,9 


a = 3, P = 2, 7 = 1 


54 


4,6 


Table 4: Transit time and absorbed dry nutrients at the end of t 


le small intestine depend on the 



values of a, j3 and 7. 



transit time was observed for a < (3 < 7, it shows the lowest level of dry absorbed nutrients. This 
stresses the key-role of "available water" W on digestion. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 1 2 3 4 5 6 



Time (hour) Time (hour) 

(c)a = /3 = 7 = 3 (d) a = /? = 7 = 4 

Figure 6: The evolution of the equilibrium Ag m — A ns depends strongly on the quantity of a, (3 
and 7. 

A s — A ns equilibrium is observed in two cases : 

(i) Firstly, it was assumed that a > (3 > 7. There is a little production of A ns from A s 
(desoubilization) because of the lack of water at the entry of the small intestine followed by 
the degradation of A s . 

(ii) Secondly, it was assumed that a < (3 < 7. There is a rapid solubilization of A ns followed 
by degradation of A s into Bi nt . The released water by degradation of A s is not sufficient 
to maintain the solution state of Bi nt since a < f3. The required water is provided by the 
transformation of A s into A ns as it can be seen in the Figure 7. 

6.3 The variation of the water ratio inside the bolus 

We are interested by the change in the ratio of W inside the bolus and its influence on the 
absorption in the small intestine. For the sake of simplicity, it was assumed that Wf ee d = W sec = 
0. 

In order to study this effect, the value of dry matter in the bolus is assumed fixed at 42 g while 
the value of water increases in each experience. 

In the numerical simulations, the ratio of water included in the bolus represented 55%, 60% 
and 66% of bolus. The ratio of absorbed dry nutrients to the total absorption (water and nutri- 
ents) is collected. Results are presented in Table 5. The numerical results of Table 5 shows that 
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(a) a > (3 > 7 (b) a < j3 < 7 



Figure 7: The equilibrium Af 71 — A ns for different quantity of a, (3 and 7. 



(%) Amount of 


(%) Ratio of 


(g) Mass of 


Transit time (hour) 


water at x = 


absorbed dry 


Absorbed dry 




inside bolus 


nutrients at x — 17 


nutrients 




50 


50 


25 


5,7 


60 


35 


26 


5,2 


66 


28 


27 


4,5 



Table 5: Dependence of absorbed dry nutrients to the water ratio. DM stands for the amount of 
dry matter in the bolus at the entry of bolus. 



increasing the value of water in the bolus decreases the ratio of absorbed dry nutrient to the total 
absorbed matter (waters dry nutrients), even though the value of absorbed dry nutrients increases. 

In fact, in our model, at the end of each experience the value of W% is approximatively 55% 
of the total mass of bolus and this equilibrium is achieved almost quickly (because of the choice 
of k w ). Therefore, increasing the value of water in this model does not have a meaningful effect 
on the final absorption and the slight increase in the absorbed mass of dry nutrients is due to 
the change of the volume of bolus in each experience which promotes the access of nutrients to 
intestinal wall for absorption, although this increase dilutes the bolus and decreases the volumic 
degradation. Here again, this is the direct result the choice of the parameters (rate of surfacic 
{ksurf) and volumic (k vo i) degradation ). 



7 Discussion 

In this paper, we have continued the modelling of the digestion in the small intestine started 
in [1]. The objective was to obtain a more realistic model of digestion process by including new 
phenomena and completing the others: DM and water are treated separately, water evolution is 
influenced by the degradation of other molecules, the effects of DF on the digestion are taken 
into account which is also a first step toward a non-homogeneous model with different types of 
feedstuffs. 

The advantages and the defects of this model as well as the perspectives are outlined in the 
following paragraphs. 

One of the main aspect of this model remains its genericity, we have tried to identify and 
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model the main generic phenomena of digestion and ignored or implicitly taken into account in 
the parameters the ones which required the specific properties of feedstuff s' molecules (effects 
on gastric emptying, on viscosity,. ..etc.). The different steps of digestion (equilibrium between 
A ns and A s , successive transformation of A s into Bi n t and B a t, s ), the effects of physical charac- 
teristic of bolus (surface and volume) on its degradation and transit, the interaction between DF 
and feedstuffs molecules (among others) have been considered while some other phenomena like 
the separation between the enzyme activities of different feedstuffs' molecules, dependence of the 
enzyme activity with respect to the dilution of bolus, the different substrates density, the impact of 
DF on initial condition (and others) have been ignored. 

To our point of view, we obtained a more realistic model by integrating these new phenomena 
in the model of digestion. In particular, the WHC of dietary fibre which in turn interferes the 
digestion of other feedstuffs molecules, leads us to introduce the separation between DM and 
water and consider all the effects of water. 

However, some other phenomena are still ignored either because of the lack of information 
concerning their effects or because of their supposed little impact on digestion at this scale. Of 
course, It would be interesting to include in the model the phenomena like interaction between 
different categories of feedstuffs and then define the specific enzyme activity for different cases to 
exploit their potential to impact the digestion. The future development of the model will be based 
on these new objectives. 

Modelling the influence of soluble and insoluble DF on the initial condition and on degradation 
of other feedstuffs' molecules as well as on the movement of bolus (Experience 6.1) is the first 
try to model a more realistic non-homogeneous bolus. DF have normally a high WHC which 
increases the volume of aqueous phase in the bolus and therefore dilutes the solution of nutrients 
and enzymes [14]. This is known to influence the volumic and surfacic reactions. However, these 
effects depend highly to the choice of W so i and Wi nso i. Another aspect of the model is the water 
equilibrium and its impact on absorption, here again the choice of parameters (W as , Wmt and 
W a bs) plays significant role on the final results of digestion. 

On the other hand, Experience 6.3 reveals the role of other parameters of model (k w , k equ i, 
ksurf, k vo i). As it has been described in this experience, the change of water absorption rate k w , 
can change ( even inverse) the numerical results. However, it is worth pointing that the choice of 
the model parameters are based on the observed behaviours (literature), the results are therefore 
consistent qualitatively with the reality (positive influence of insoluble DF on digestion, negative 
effect of soluble DF, ...). 

Taking into account these new phenomena requires the introduction of new parameters which 
can be identified with the help of existing experimental data. To the best of our knowledge, some 
of these parameters like the water associated to dry nutrients (W so i , Wi nso i, W as , Wi n t and W a b s ) 
are introduced for the first time and they should be identified 

One of the advantages of this model is its capacity to be reduced and to be adapt to the new 
experiences which makes the parameters identification possible. Reduction consists in the decreas- 
ing the number of equations of system or the number of parameters without affecting its genericity 
(e.g. a bolus which does not contain the DF, results in a more simplified digestion process which 
in turn caused decreasing the equations (parameters) involved in the digestion model). 

It is also worth pointing that the value of most of the parameters depends to the special category 
of feedstuffs. The close collaborations between biologists and mathematicians is therefore required 
to identify these new parameters (literature data in biology, define the new experimentations, ...). 
This reveals one of the main interest of modelling which is to ask the precise questions about the 
modeled phenomenon. In fact, this approach allows to use all the existing data and limit the new 
animal experimentation to the special cases (when the existing data are not sufficient). 
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Appendix 
1 The Model equations 

o Equation of transport 

d 2 x _ j dx cp + c\r so i K visco dx 

dt 2 ~ T{ dt> a + bx [W] dt' 

o Non solubilized substrate A ns 

A A 

-^ = -k equi {fi([W])A ns -A s ) 

o Soluble DF F so i 



dFt 



dt 

o Solubilized substrate A d s m 

j A dm 

an. a , , w dm 



so1 - -k s e exoP h(x)[F^]V app 



dt 



k vo i{x)[Af n ]V app - k surf [A s ][W]S so i + k equi (n([W])A ns - A s ) + secretions 



o Intermediate substrate Bf^J: 

=k vo i{t)[Af n ]V app + k s e exo ph{x)[F^]V app - k sur f[B int ][W]S so i + secretions 
o Absorbable nutrients 

" " hs = k surf [B int ] + k surf [A s ] [W]S sol - k abs [BiZ\S S oi 



dt 

o Exogenous enzymes e(t) 



o Free Water Evolution W(t) 

dW 



de 
dt 



dt 



ak equi {n{[W])A ns - A s ) + (a - (3)k vol (t)[A<!: m }V 



app 



surf [Bint] + (a - 7)k surf [A s 
+ (A, - mse eX oph(x)[F^]V app + jk abs [B±]S sol - k w (W - QAM). 



o Evolution of the bolus mass M(t) 



= WZw^- kw ^ W ~ °- 1M) " k abs[Baf s ]S sol + Secretions) 



o Evolution of bolus volume V(t) 

dV 



= -k w (W - 0AM) + Secretions 
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